Impact of Quantum Phase Transitions on Excited Level Dynamics 
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The influence of quantum phase transitions on the evolution of excited levels in the critical pa- 
rameter region is discussed. The analysis is performed for ID and 2D systems with first- and 
second-order ground-state transitions. Examples include the cusp and nuclear collective Hamiltoni- 
ans. Applications in systems of higher dimensions are possible. 
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I. INTRODUCTION 



In the last decade, the concept of quantum phase tran- 
sitions (QPT) 0, @, Q has triggered a lot of activity 
in solid-state physics, see e.g. Refs. 0, la, ]g|, and in 
many-body (nuclear) physics, see e.g. Ref. [7| and the 
references therein. The QPT shows up in the system's 
infinite-size limit at zero temperature as a nonanalytic 
change of the ground-state properties with an external 
control parameter (e.g. an interaction strength). Exam- 
ples range from order-disorder transitions in spin lattice 
systems [4| to transitions between various quasidynami- 
cal symmetries in some interacting boson/fermion models 

Phase transitional effects were recently identified also 
for excited states in the two-level pairing models that ex- 
hibit a second-order ground-state QPT |H Q E [3 . 
Two characteristic signatures of excited-state transitions 
were recognized, namely (a) an anomalous evolution of 
individual excited levels as they cross the phase separa- 
trix in the plane of the control parameter versus energy, 
and (b) a singular behavior of the level density at the 
critical energy [lq ]. In this article, a generalization of 
these results to systems with the first-order ground-state 
QPT is discussed. It is shown that a phase transition 
at zero temperature has significant consequences for the 
dynamics and density of excited levels in the critical re- 
gion of the control parameter. These features are closely 
related to thermodynamical properties of the system. 

The plan of the paper is the following: In Sec. |TT] 
we introduce a general framework for studying excited- 
state phase transitions in quantum systems with coincid- 
ing classical and thermodynamical limits. This type of 
systems underlies a wide class of models used in many- 
body physics. In Sec. IIII1 the concept of excited-state 
QPT is illustrated by a simple, prototypal example of 
a one-dimensional model with both first- and second- 
order ground-state transitions (the cusp catastrophe). 
The dependence of phase-transitional effects on the di- 
mension is discussed in Sec. IIV1 where an example of 
a two-dimensional system, closely related to the nuclear 
collective model, is analyzed. The last section contains a 
brief summary and conclusions. 



II. QUANTUM PHASE TRANSITIONS FOR 
EXCITED STATES 



As mentioned above, quantum phase transitions are 
related to the zero-temperature limit of the system, thus 
to properties of the ground state. An extension of the 
QPT analysis to excited states can be achieved through 
the free energy 



F = Tr(gH) +T Tr(g In g) = -TlnZ, 

(E) ~~^S~~ 



(1) 



where H is the Hamiltonian, g the canonical density op- 
erator, Z the partition function, (E) an average energy 
at temperature T, and S the corresponding entropy [17j |. 
For T — 0, the free energy coincides with the ground- 
state energy. 

The known relations 
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where ((E 2 )) stands for a thermal dispersion of en- 
ergy, set up the conditions for thermal phase transitions: 
While in the first-order transition the entropy suddenly 
jumps at a certain "critical" temperature T c , implying a 
discontinuous first derivative of F, in a continuous phase 
transition only the second and/or higher derivatives are 
affected. Such situations can only occur in the limit of 
infinite system size, when the canonical description is as- 
sumed to converge to the microcanonical one. Then the 
entropy can be written as S ex lnp, where p is the level 
density at energy E = (E). 

If the Hamiltonian depends on an external control pa- 
rameter A, so does the free energy ((T|). An interesting 
question is whether a quantum phase transition at T = 
and A = A c (0) extends to T > in the above thermo- 
dynamical sense, forming a phase separatrix A C (T). Al- 
though specific examples exist, in which the T — critical 
point is isolated, in generic situations the phase transition 
exists in the XxT plane [H S S 0> @] • The classification 
of such transitions is again in terms of the behavior of the 
free energy at the transitional point, but with respect to 
both variables T and A. In particular, the first-order or 
second-order transitions, respectively, imply a disconti- 



nuity of the gradient vector or of the curvature matrix 
(Hessian) associated with F(X,T). 
Assume a linear dependence of the Hamiltonian 



H(X) = H + XH X 



(3) 



on a single real parameter A, with Hq and H± standing 
for two incompatible dynamical terms. From the pertur- 
bation theory it follows that the additional relations to 
@ read as: 
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Here ((XY)) = (XY) - (X)(Y), with (•) standing for 
the thermal average of individual level energies Ei and 
their derivatives E t = j^Ei, Ei = 4^Ei. 

The formulas in Eq. @ can be used to anticipate the 
implications of thermal phase transitions on the level dy- 
namics. The first-order transition, which shows up as a 
discontinuity of -§^F and -§\F, leads to a jump of the av- 
erage slope (E) of energy levels. If the phase separatrix 
is not parallel with A or T, the jump of (E) is connected 
with a jump of p, thus also a jump of the microcanonical 
entropy S. This may be viewed from an analogy with the 
ray refraction on a tilted interface — the vertical distance 
of rays (alias level spacing p~ x ) changes at the interface. 
On the other hand, the second-order transition, with dis- 

r) 2 r) 2 

continuous -SfjF and -§vzF, is linked to changes of the 
average level curvature (E) and/or of the slope dispersion 
{(E 2 )}. Singular derivatives generate continuous phase 
transitions with no Ehrenfest classification. 

In the following, we will illustrate the above general 
findings by concrete examples. We focus on systems with 
a finite number of quantum degrees of freedom /. In the 
many-body context, such systems usually arise from bo- 
son/fcrmion models based on the dynamical Lie algebras 
of finite dimensions [18 ] . A typical example is a system 
consisting of a conserved number N of interacting bosons 
with the single-particle Hilbert space of dimension n: in 
this case / = n — 1, with U(n) being the dynamical al- 
gebra [ig |. 

A significant feature of the above type of systems is the 
fact that the infinite-size limit coincides with the clas- 
sical limit @, y, 0, |2Cj . In bosonic systems, e.g., the 
value of A -1 can be associated with the Planck constant 
ft,. Therefore, the asymptotic-size limit of the level den- 
sity p(E), which constitutes various kinds of excited-state 
phase transitions, is proportional to the classical phase- 
space volume V{E) available at energy E, 

V{E)= 1 5(E~H) d f p d f x = — 6(E - H) d f p d f x . 

W(E) 

(5) 
Here, H stands for the classical Hamiltonian, 9 is the step 
function, and W(E) represents the phase space volume 
available at energies less than (or equal to) the value E. 



This formula will be employed to classify the excited- 
state QPTs in the forthcoming examples. 



III. CUSP HAMILTONIAN 

Consider first a one-dimensional Schrodinger equation 

(6) 



2 dx 2 



H = — -y-j + x* + ax 2 + bx , 



with a potential having the well known cusp form [2l|. 
Here, a and b are the cusp parameters and K = -4= is 
a classicality constant bearing information on the Planck 
constant and on the mass parameter M. As discussed 
above, in the bosonic models h is inversely proportional 
to the size of the system, measured by the boson number 
A, the effective mass depending on a concrete applica- 
tion. Note that the Hamiltonian |6]) can be transformed 
[2(j | to a form equivalent to the model containing two 
types of interacting (pseudo)scalar bosons (the bosonic 
formulation of the Lipkin model [22j,[23j]), although in the 
latter case the kinetic part is generally more complicated 
than that considered here [la]. While the difference in 
kinetic terms obscures a direct comparison of K and N, 
it is clear than the limit K — > corresponds to N — > oo. 

The cusp potential (with germ a; 4 ) represents the most 
common type of catastrophe [211 ] in dimension one which 
generates both first-order and second-order phase tran- 
sitions. These can be associated with two trajectories in 
the plane a x b, namely with potentials V\ — x A — x 2 + Xx 
and Vi — x A + Xx 2 depending on a single parameter 
A. The potential V\ has two minima at x =/= within 
the region demarcated by a pair of spinodal points at 
A = ±— -=. If A in V\ varies from negative to posi- 
tive values, the ground state exhibits a first-order phase 
transition — the swap of both minima at the "critical" 
point A c = 0. On the other hand, if A in V2 varies 
from positive to negative values, the ground state shows 
a second-order phase transition — the single minimum at 
x = (A > 0) splits at A c = into a pair of degenerate 
x =/= minima (A < 0). 

It is clear that the cusp potential generates the above 
types of ground-state QPTs only in the semi-classical 
limit K — > (equivalent to N — > 00 in the bosonic formu- 
lation), when the zero point motion vanishes. The ques- 
tion is what happens in this limit with excited states close 
to the critical point. Let us first analyze the behavior of 
the phase-space volume ([5]). It turns out that V(E) has 
some non-analytic features, which are connected with the 
bottom energies E^ and E^> of both potential wells and 
with the top energy E 1 - 3 ' of the barrier in between. For 
a particle moving in a parabolic minimum with classical 
frequency u> we find V = — . Therefore, for the double 
well potential (both minima are locally quadratic) the 
phase-space volume V(E) starts from a non-zero value at 
E = E^ 1 ' and jumps to a higher value at the energy of 
the secondary minimum E^ 2 ', where new states become 




FIG. 1: Level dynamics in (a) first-order and (b) second-order 
phase transition generated by the cusp Hamiltonian ([6]) with 
(a, b) — ( — 1, A) and (A,0), respectively; K = 10 -2 . The inset 
in panel (b) shows the subset of levels with positive parity. 



available in the second well. This results in a first-order 
phase transition. 

The excited level dynamics for the cusp potential V\ 
is shown in panel (a) of Fig.[TJ Note that the numerical 
diagonalization was performed in a truncated oscillator 
basis, the convergence issues being fully under control. 
The V-shaped region, which is apparent in the figure, 
coincides with the A x E domain where states can be 
located in both potential wells. As levels hit this domain, 
their "laminar" flow turns into a "turbulent" one. Note 
that although the mutual approach of levels in the V- 
region is rather sharp for low-lying states, the crossings 
are always avoided due to the tunneling effects (see the 
inset). As the slopes of up- and down-going stretches 
of an individual level trajectory compensate each other, 
the average slope changes abruptly in transition to the 
V-region. Indeed, this leads to a sudden increase of the 
level density, as anticipated in the previous discussion. 

The zig-zag pattern of level trajectories developed in 
the V-region, Fig.QJa), is connected with two families of 
states localized in the first and second potential well. Un- 
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2: Squared wave function of the 5th state from Fig.QTa). 
\ip axis is not drawn. 



der the neglect of tunnelling, these states form approx- 
imate eigenstates of the Hamiltonian. While states in 
the up-going well rise in energy, those in the down-going 
well decline. As an eigenstate propagates through this 
region, its wave function alternates between the two lo- 
calizations, with the changes taking place at each avoided 
crossing. This is illustrated in Fig. [2 where one observes 
that the structures with 1, 2, 3 etc. peaks cross the gaps 
between neighboring levels and continue along the direc- 
tion in which the respective well moves. Interestingly, 
the wiggling pattern of energy levels inside the V-region 
becomes infinitely dense (undifferentiable) in the semi- 
classical limit. 

On the top E — E^ of the quadratic barrier separat- 
ing both minima of the cusp potential the phase-space 
volume V(E) exhibits a logarithmic singularity. It is 
connected with asymptotic dwell times of classical tra- 
jectories at the maximum which yield a locally vanishing 
gap between quantum levels. This singularity was ana- 
lyzed e.g. in Ref. [2J]. Its impact on the level dynamics 
is observed at the upper side of the V-region in Fig.QJa). 
When level i reaches the energy E^\ the curvature Ei 
becomes discontinuous and infinite for K — > 0. This leads 
to a locally infinite growth of the level density and implies 
a kind of continuous phase transition (with no order). Let 
us note that the decrease of p(E) after the singularity is 
connected with a concave increase of W(E), see Eq. (0, 
which indicates anomalous thermodynamical properties 
[17| . As shown below, this peculiarity is specific for ID 
(or quasi-lD) cases. 

Panel (b) of Fig.Q] shows the level dynamics for the 
second-order phase-transitional potential Vi . The behav- 
ior observed is essentially the same as that in the Lipkin 
model without parity violating interactions. All levels 
are characterized by the parity quantum number and for 
£<0 they form approximately degenerate parity dou- 
blets. The phase-transitional evolution is detected at en- 
ergy E = E 1 - 3 ' — corresponding to the local maximum 
of the potential at x = 0. This transition has been stud- 
ied in the Lipkin model [25[ as well as in a wider class of 



two-level pairing models used in nuclear and molecular 
physics [3, [ljj [la, [26| • Let us note that the thermody- 
namical description of this effect, based on the mean-field 
approach, was presented in Ref. |2j. 



IV. 2D COLLECTIVE HAMILTONIAN 

The cusp catastrophe applies to ground-state QPTs in 
a large class of models, including e.g. two-level inter- 
acting boson models [20J. However, since the number of 
degrees of freedom / is usually larger than 1, the excited 
spectrum differs from that in the cusp example. 

For / > 1, only a fraction of trajectories may cross the 
stationary point, hence a general trend is a smoothening 
of non-analytic features at the critical energy. In partic- 
ular, for / = 2 only -j^V evolves in a discontinuous way 
when crossing the energy of a local minimum or maxi- 
mum of the potential, i.e. the transition is of the second 
order in both cases. For a saddle point, the derivative 
^5 V is "continuous" but infinite (V has a singular tan- 
gent). A further increase of / leads to even higher orders 
of transitions. These statements can be verified by ex- 
plicit evaluation of Eq. ([5]) for the respective forms of the 
potential in general dimension. As a consequence, for 
/ > 1 one obtains continuous (second-order or softer) 
phase transitions along all three sides of the V-shaped 
region of phase coexistence in the first-order QPT. 

To illustrate these conclusions, we analyze the follow- 
ing / = 2 Hamiltonian 



H 



2 



ld_ d_ 

r dr dr 



i d 2 



dip 2 



+ r 4 + ,4r 2 +Br 3 cos3</?, 

(7) 

where (r, <p) are polar coordinates and {K, A, B} ad- 
justable parameters (K has the same meaning as in the 
cusp case). We consider only the states with a periodic 
boundary condition on the sextant ip e [0, ?]. This sys- 
tem is closely related to the geometric model and the 
interac ting boson model of collective motions in atomic 
nuclei [ijj, [27( . These models describe 5 degrees of free- 
dom related to nuclear quadrupole deformations, i.e. 2 
deformation parameters (/3, 7) = (r, p>) and 3 Euler an- 
gles, but for states with zero spin the Euler angles be- 
come frozen. The potential reads as that in Eq. |(7J). Let 
us note the genuine 5D nuclear Hamiltonian for zero-spin 
states has a slightly different kinetic term than the 2D 
Hamiltonian (|7|), but both excitation spectra show the 
same qualitative features [281 ]. 

For K — > 0, the Hamiltonian ([7]) exhibits a QPT along 
the parabola B 2 = AA that separates two basic forms of 
the potential. The transition is of the first order except 
at B = where it is of second order. On the inner side of 
the parabola, the potential describes a quartic oscillator 
with the global E = minimum at r = 0. On the outer 
side, the potential with ip e [0, %) has the global E < 
minimum at r > 0, <p = ? (for B > 0) and a saddle point 
(with E < 0) at r > 0, p — 0. In the spinodal region 
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FIG. 3: Level dynamics for Hamiltonian (J7J) close to the first- 
order phase transition (A = X, B = 1, K — 10 -3 ). Dashed 
lines demarcate the phase-coexistence domain. 



of the first-order transition, the r > minimum coexists 
with the r = one, creating another saddle point (with 
E > 0) in between (at ip = §). For A < 0, the E = 
minimum at r = turns into a local maximum. 



We consider a scaled potential V 



-Xr 2 



' cos 3ip 



with the critical point A c = 4 and two spinodal points 
at A = and Jj. The level dynamics in the phase- 
transitional region, obtained by a numerical diagonal- 
ization of the Hamiltonian in a modified 2D oscillator 
basis, is depicted in Fig.[3l Computational details will 
be given elsewhere 28] . As in Fig. [Ha) , all levels pass 
through virtual (avoided) crossings. The legs of the 
phase-coexistence "triangle" (dashed lines) correspond 
to energy E^ 2 ' of the secondary potential minimum (a 
second-order phase transition) and the upper side E^> 
represents a saddle point (a continuous phase transition 
with a singular tangent of V) . 

Despite some similarities, the picture is less dramatic 
than that for the ID cusp Hamiltonian. Under a scrutiny 
one may observe that the refraction of levels along the 
boundaries changes the dispersion of slopes rather than 
their average, consistently with the softer types of phase 
transitions. Note that for B = one can separate sub- 
sets of levels, characterized by the 2D angular momen- 
tum quantum number m, with gradually lowering phase- 
transitional signatures [lj, LL6| , but this cannot be done 
in general as m is not conserved for S^0. 

We have evaluated the level density from the numeri- 
cal spectra at selected values of A and show the result in 
Fig. [4] Parameter K, which scales the absolute density 
of states, was chosen for each A to ensure a sufficiently 
dense spectrum. While the non-degenerate double-well 
potentials (A = 0.22, 0.27) generate typical forms with 
three characteristic slopes of p in the regions E < E^ 2 \ 
E^ < E < E<< 3 \ E^ < E (see the marks), the critical 
potential (A = 0.25) yields a shape with only one transi- 
tional point. Changes of p at the critical energies become 
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FIG. 4: The level density for Hamiltonian Q with (A, B) - 
(A, 1). Parameter K was set separately for each histogram 
within the bounds 10" 4 < K < 5 ■ 10" 3 . The binning of 
histograms is shown by small dots. Marks indicate minima 
(T), maxima (a) and saddle points (■) of the corresponding 
potentials. 



sharp in the K — ► limit. The level density above the 
coexistence region (A = 0.5) exhibits no transitional fea- 
tures. In contrast, the spectrum below the lower spinodal 
point (see the inset) displays transitions connected with 
the (p = saddle point and with the r = local maxi- 
mum. 

Apart from indicating the phase structure of the spe- 
cific model under study, the curves in Fig. 3] illustrate the 
dctectability limits for continuous phase transitions in a 
finite system. Let us stress that the transitions would be 
even more smoothened in higher dimensions. An appli- 
cation of these results in the interacting boson model (for 
the zero-spin spectrum) is in order. 



suits for first-order ground-state transitions which sup- 
plement earlier findings for the second-order transitions 
[lj, LL5J, LL6|. In dimension one, the cusp Hamiltonian ((6|) 
is a fundamental example of both first-order and con- 
tinuous transitions. Despite the fact that the relevant 
physics in this case is just the basic-level quantum me- 
chanics, the effects observed constitute the clearest real- 
ization of quantum phase transitions affecting individual 
excited states. 

For higher dimensions, the signatures of excited-state 
phase transitions are weakened in a twofold sense: (a) 
The transitions only affect higher derivatives of the level 
density, and (b) only some bulk properties of the level 
dynamics are influenced rather than evolutions of all in- 
dividual states. These features — which on the classical 
level can be linked to an increasing size of the phase 
space, hence a decreasing impact of singular (e.g. station- 
ary) points on the classical motions — hinder the practical 
detection of excited-state QPTs in finite samples of the 
system. 

Finally, let us stress that the findings discussed in this 
article are relevant for the systems with synonymous 
infinite-size and classical limits. This feature allows one 
to associate the asymptotic level density with the phase- 
space volume and therefore to unambiguously classify 
the excited-state phase transitions. Numerous examples 
of this type of systems can be found in nuclear, molecu- 
lar, and mesoscopic physics. It would be interesting to 
learn how the features discussed here extend to the infi- 
nite lattice systems that do not share the above property. 



CONCLUSIONS 
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